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CHARACTERISTICS OF TURBULENCE IN A BOUNDARY LAYER 
WITH ZERO PRESSURE GRADIENT ‘ 

By B. S. Kj.kbanoff 


SUMMARY 

The results of an experimental investigaiUm of a turbulent 
boundary layer mith zero pressure gradient are presented. 
Measurements with the hot-wire anemometer were made of 
turbulent energy and turbulent shear stress, probability density 
and flattening factor of u-fluciuation {fluctuation in x-direction) , 
spectra of turbulent energy and shear stress, and turbulent 
dissipation. The importance of the region near the wall and 
the inadequacy of the concept of local isotropy are demonstrated. 
Attention is given to the energy balance and the internnftenl 
character of the outer region of the boundary layer. Mso 
several interesting features of the spectral distribufion of the 
turbulent motions are discussed. 

INTRODUCTION 

Tlie statistical theory of turbulence introduciul by Taylor 
and elaborated upon by Karinan, Howarth, and Dryderi lias 
played an important role in provdding a sound basis for tlu^ 
study of turbulence. In the main, this advance lias been 
confined to the domain of homogeneous and isotropic 
turbulence. At present turbulent shear flows present difli- 
culties so formidable that statisti(^al theories have made little 
progress. Tho older semiempirical and phenomenological 
theories still constitute the most tangible theoretical met hods. 
It is now generally believed that experiment slioiild lie called 
upon whenever possible to furnish information on the actual 
behavior on the grounds that such information is needed to 
acquire an insight into tlie turbulence processes and to form 
the basis for a sound tiicoretical approach. 

Experimental investigations in shear flow have made 
encouraging progress, largely because of the increasing num- 
ber of statistical properties that may be measured by im- 
proved hot-wire techniques. Measurements have been 
made in the jet, wake, two-dimensional chanmfl, pipe, and 
boundary layer (e. g., refs. 1 to 8). These have shown that 
the basic assumptions of the phenomenological theories are 
inconsistent with the experimental evidence as to tlie nature 
of the turbulent motions. The work of Corrsin (ref. 1) 
revealed the intermittent character of turbulent flow near a 
free boundary. This was later studied in sonu^ detail by 
Townsend (ref. 2) and is now recognized as a phenomenon 
associated with a sharp but irregular division between tur- 
bulent and nonturbulent flow. The recent theoretical 
contributions of Kolmogoroff (refs. 9 and 10), Heisenberg 
(ref. 11), and others dealing with the concept of local iso- 
tropy and the spectrum of turbulent energy have givim added 


impetus to experimental work in sliear turbulence and have 
encouraged the point of view that some of tlie properties of 
isotropi(‘ tuihuhuH*!' may Ix* appli(‘ab]e. Tlie application 
of isotropy to sluxir flow has been ratlier extensively studied 
in wakes (nd. 8) and found to be useful, but its usefulness in 
a boundary layer is still cj[uestionable, especially in tlie 
region near tlu^ wall. At best, isotropy belongs to the final 
or dissijiation slagi* of shear turbulence and here can throw 
little light on the mechanisms pertaining to sliear flow. Any 
attempt to invi^stigate the turbulence mechanism is ham- 
pered by the lack of an exjierimental technique for measure- 
ment of the ])ressure flinduations. However, it is possible 
to obtain information pertaining to (he energy balance, the 
study of which may Ix^ considered a propiu* approach to the 
problem. An inijiortant step in this direction in the investi- 
gation of the boundaiy layer was made by Townsend (ref. G). 
With respect to this jioint of view tlie jiresent work attempts 
to ])rovide additional information, especially for the region 
close to the wall, and to obtain a more direct measure of the 
dissipation without relying comjiletely on the concept of local 
isotro])} . 

Attention her(‘ is given to the boundary layer with zero 
pressuiT' gradient undiu* conditions as favorable as practicable 
for the use of hot-win* technique, ddie method of obtaining 
tin* boundary layer is described in reference 8, and the present 
work may be i*(‘gai(h*d as a continuation of the former work 
in whi(*h use* is now made of an artificially thickened layer as 
a research tool. While there are many important features 
of boundary-layer turbulence, the present investigation is 
concerned with three major phases: The intermittency and 
its effects at the free boundaiy, energy balance, and the 
spectral distribution of turbulent energy and shear stress. 
Since the boundaiy layer is a complex flow combining the 
effects of a free boundary on the one side and those of a solid 
wall on the other, no one part of the layer could be de- 
(‘Uiphasized. This makes it an interesting subject if not a 
simple one. 

The present investigation was conducted at the National 
Bureau of Standards under the sponsorship and with the 
financial assistance, of the National Advisory Committee for 
Aeronautics. The author wishes to express his appreciation 
to Dr. G. B. Schuhauer for his active sujiport and construc- 
tivi! criticism and to ai'knowledge gratefully the assistance 
of Miss Z. W . Diehl and Mr. T. J. Kelly. He also wishes to 
thank Dr. (h AI, '^rchen and Dr. J. Laufc'r for their interest 
and stimulating discussions. 
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SYMBOLS 
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coofficient of skin friction, To/4 pV 
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K{k) three-dimensional spectral function associaU'd 

with k 

e voltage fluctuation 

signal produced hy ^/-component 
^2 signal produced by z;-coinponent 

F flattening factor __ 

percent of turbulent energy associated with k\ 

fuiki) turbulent energ}' associated with ky 

f7v{k\) turbulent shear stress w associaUal with k^ 

tui'bulent energy F associated with k^ 
k tliree-dimensional wave number 

k,=(:^xhl8yY^ 

one-dimensional wave numbtu' 
n frequency 

P{u) I)rohability density of u 

P{uiu') probability density of 7//u' 

, ,, , duv dlJ 

mean static pressure 
mean static pressure at j:=10.r) feet 
instantaneous pressure fluctuation 
reference dynamic pressure ahead of plate 
longitudinal space correlation (‘oefficient of 
^/-fluctuation 

lateral space correlation coeflicituits of i/-fluctua- 
tion 

skewness factoi* 
time 

r-component of mean v(‘locity 


P 

Vo 

Vi 

<io 


R, 


s 

t 

V 


’ A p 

?/, V, ‘W 

7/', V\ W' 

7/h r, w‘ 
1W 

V 


mean velocity in free stream 

instantaneous turbuhuit velocity fluctuations 
in x-y ?/-, and --directions, r(‘si)eclively 
i‘oot-mean-s(iuar(' valu(‘S of ?/, r, and w 
imain-square values of 7/, v, and w 
tui'bulent shearing stress 
7 /-(‘omponent of mean velocity 


x„ distance along surface from virtual origin of 

boundary la}^er 

y distance normal to surface measured from sur- 

face 

' V 

z direction perpendicular to X7/-plane 

T intermittoncy factor 

6 boundary -layer thickness 

be boundary-layer energy thickness, 

c dissipation of turbulent energy in isotropic 

tui'buhuice 
f=(V2a/5)-> fi//6)-0.781 

B angle between velocity ve(*tor and normal to 

wire 

longitudinal microscale 
V kinematic viscosity 

p density of air 

(T standard deviation 

Ty shearing stress at wall 

X ahsolute constant 

coj. ar-{‘omponent of vorticity 


( ) 


mean value 


PROCEDURE AND RESULTS 

EXPERIMENTAL ARRANtJEMENT 

The pn^sent investigation was conducted at the National 
Bureau of Standards in the 4^-foot wind tunnel shown in 
figure 1. The turbulence level of the tunnel was 0.02 per- 
cent at ^10 f(‘ct j)er second and 0.04 percent at 100 feet per 
sex'ond. The low level of tiirhulence was obtained by damp- 
ing screens placed in the settling chamber. The boundary 
layer was develop('d on a smooth, flat, aluminum plate 12 
feet long, 4 I 2 hnd wide, and % inch thick with a symmetrical 
and point ( h 1 [(aiding edge. The plate was mountcal vertically 
and centrally in the test section of the tunnel. 

The scheme of artificially thickening of tur})iilent bound- 
ary hqyers developed in nderence S was applied herein order 
to realiz(^ tlu^ advantages of a thick boundary layer, namely, 
the larger s(‘ah's of mean and fluctuating flow fiedds which 
d(M‘reas(* errors due to th(3 finite size of the hot-wire proh(' 
and the limited upper-fi'equency r(*spons<‘ of the equipment. 
In addition, high Reynolds numlxu's are afforded without 
high speeds, d'lu* tlii(*k(*ning was a(‘hiev(‘d by (*ov<u'ing tin* 
first 2 feet of the plate with sand roughness consisting of 
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No. IG floor-sanding paper. At the working station lOK 
feet from the leading edge tlie boundary layer was 3 inches 
thick. All of the measurements were made at this posi- 
tion. The frec-stream speed was 50 feet per second. The 
elaborate tests described in reference 8 showed that the 
boundary la}^er was tlie fully developed, smooth-wall type, 
having an apparent d(*velopment length of 14.2 feet of smooth 
surface. The correspoiiding lengtli Reynolds number based 
on values of x measured from the virtual origin was 4.2 X 10®. 

In order to obtain a condition of zero pressure gradient 
along the ])late, the passages between the tunnel wall and 
the platen was made sufficiently divergcuit to offset tlie na- 
tural fall in pressure du(^ to boundary-layer growth. This 
was accomplished by a flexible side wall whicli could be 
j)ositioned i)y s(*rews threaded into the tunnel wall. The 
final pressure distribution is shown in figure 2. From 2 feet 
on downstream tlie pressure is seen to be constant on the 
average with variations about ]i percent of qo about the 
mean. These variations were associated with the inlierent 
waviness of the surface and could not be removed by the 
adjustable wall, 

INSTRUMENT MOUNTING AND TRAVERSING 

The positioning and moving of an instrument from point 
to point was always done in the manner that best suited the 
instrument and the purpose. For example, the pressure 
distribution along the plate was determined by means of 
a static tulx^ mounted from a carriage that could be moved 
and j)ositioned by remote control. Since all other types of 
measunumMits w(u*e made at the 10.5-foot station, the various 
measuring probes were supported on rods extending through 
tlie plat(‘ to a micrometer-screw trav(‘rsing device on the 
opposite side. This provided adequate rigidity and negli- 
gible interference and permitted movement by known 
amounts to and from tlu* surfa(‘e. Initial distance from the 
surface was obtained by using a prism to reflect the images 
of tlu‘ suilace and the probe on the calibi-ated scale of a 
inicros(‘ope. 

In special (*as(‘s such as the measurtumuits of the Hy and 
corndations, small manually operated ti’aversing units 
that could be mounted to the rods were employed. In all 
cases adecpiate rigidity and freedom from interference rather 
than (‘onv(uiien(‘e of op(‘i*ation ditdated the arrangement. 

MEASUREMENT OV PRESSURE AND MEAN VELOCITY 

Nickel tubing 0.04 inch in diameter and 0.003 inch in wall 
thickness was used for the impact and static-pressure tubes. 
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Figckk 2. — Pr(‘ssuro clisiribul ion along llal plato. 


The static-pressure tube was made according to the coir 
ventional design for such a tube, and the pressure distribu- 
tion was measured by traversing longitudinally at a dis- 
tance }i inch from tlie surface. The impact tube was flat- 
tened at the end to form a rectangular opening 0.014 inch 
wide and, together with a static tube similar to that used 
for the pressure distribution, was used to measure the mean- 
velocity distribution. Velocity distributions were also 
measured at 10 inches above and below the center line. 
These agreed well with the distribution at the center line 
and thus confirmed the two-dimensional nature of the flow. 
No correction was made for the effect of turbulence on the 
measured values nor for the effective geometric center of the 
impact tube. 

Close to the wall, 0<^?/<C0.05 inch, the mean velocity was 
measured with a hot-wire. Platinum wire 0.0003 inch in 
diameter and approximately % inch long was used and 
operated at low-temp(U’ature loadings in oi*der to minimize 
the influence of the wall on the heat-loss characteristics of 
the wire. A correction was made for the effect of the turbu- 
lence level on the measured hot-wire values. This was done 
by a graphical method using the mean-velocity-voltage cali- 
bration curve and the measured root mean square of the 
voltage fluctuations. The corrected mean velocity was 
higher than the observed, witii a maximum correction of 
about 10 percent. The hot-wire values were in good agree- 
ment with the pitot-static-tube values. 

The velocity distribution is shown in figure 3. The dashed 
line denoting the velocity gradient at the wall was computed 
from the shear at the wall and is in satisfactory agreement 
with tlie measured values. 

TURBULENCE INTENSITIES AND SHEARING STRESS 

The hot-wire equipment used is described in detail in 
reference 12. The frequency response of the uncompensated 
amplifier was flat from 2 to 50,000 cycles and was down 17 
percent at 70,000 cych^s. A number of cutoff filtius were 
provided to cut out frequencies above tbe range needed in a 
particular measurement so that unnecessary noise was 
eliminated. By this expedient the input noise level was held 



Fi<;irKP: 3. — Moan volooity distritnition. 
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down to tlio order of 2 to 4 microvolts, (i'ompensation for 
the time lag of tln^ wire was determiruHl by the square-wave 
method and was accomplislu'd by a resistan(*e-eapacitance 
network in the amplifier. Platinum and platinum-rhodium 
(90 percent plaliriiim and 10 percent rhodium) wires 0.0001 
inch in diameter and % millimeter long having n time constant 
of approximately 0.25 millisecond were generally used. In 
cases where the signal-to-noise level pn‘sented some* difficulty, 
0.00005-inch-diameter wire was used. This was wire drawn 
by the AVollaston process and the method of attaching iIh' 
wire was to etch away the sihuu* and to soft-solder the wir<‘ 
to the supporting ])rongs. Fim‘ sewing ti<hm1I(^s wimh' us(‘d as 
the prongs for the u-hohha* and fine jewc'hu-’s broaches wei't' 
used as prongs for the X-wii-es. No wire-length corivctions 
were applied to anv of the data. T\w method of measuring 
the turbulence intensities ?//, v\ and w' is describe<] in reh'r- 
enc(' Ph u' , v', and ?// (list I'ibut ions ai'e shown in lignin' 

4. The value of u/ was obtaimul as close to the wall as 0.004 
inch, but because of the com])aratively larger size of the 
probes necessary for the measui-enK'ut of r' and w/ it was 
not possible to measui*e tlu'se closer than al>out 0.045 inch. 
The exti*a])olated values of r' and w' wen' obtained by (‘om- 
paring the boundary-layer valiu's with those obtaiiH'd in a 
pipe (i*ef. 5) on a basis of Ly against //*. Good agn'ement 
existed within the measuiTd nuige and cons(*quently the pip(‘ 
data which contained valiu's coiTcsponding to distanc'cs closer 
to the wall (because of the lower* Reynolds number) wc're used 
as a basis for exti*apolation. 

The measiiiTinent of iTr was made by tlu' conventional 
X-wire method (h'serribed in i*eference 15. Use was inadi' of 
the experimental data for the angle rr'sponse found by 
N('wmari and Leary (ref. 14), nanudy, (cos rather than 
the customary (cos wher-e 6 is the angle l)etween th(' 
velocity vector and the normal to tlu^ win'. The signals to 
be dealt with are Cj and C 2 produced by tlu' u- and iM-om- 
ponents, i-espectiveh'. In princi[)le, oru' wire of an X-wire 
pi'obe contr'ibutes the output {Vi^C 2 )“ and tin' other* con- 
tributi's {ex — e^y^. The difference l>etween tlu'se two outputs 



giv(‘s €^€.2 fi'orri whic'h uv may be calculated by empkyving 
th(' expc'rimentally detei’iniried faxdors of proportionality. 

The rc'sults arc shown in figure 5. Since' the viscous shear- 
ing stre'ss reache'd oidy 2 percent of the total at the point 
ne'aiest the wall, the turbulent shearing stress uv shows the 
characUi'ristics of the total. It appi-oaches the wall with 
z('ro slope as it must when the pi*essur(‘ gradient is zero and 
is in ('xcellent agi‘(*ement with th<‘ value at the wall cah'ulatc'd 
by the method of Squir*e and Young (r*ef. 15). 

SI’ECTKA OF u2, v^, AND 

Tile ni('an-s(piare values of the u-flucluation Ix'tween IVi'- 
(pu'ncies }\ and n br/;? were obtained by fi'cding the signal 
from tlu' (‘ompi'usated amplifier into a General Radio wave 
analyzer having a frequeruy range from 10 to 10,000 cps and 
a fix('d band pass. Offie signal from the arialyzer was then 
f('d into a t lu'rmocouple circuit for* measuring tlu' int'an- 
square output. In reading the output, averages w('re tak(*n 
for a period of 1 minute. 44 h' mi'an-square voltage asso- 
ciated with th(' hot-wire signal passed by th<‘ band was ob- 
taiiu'd by h'i'ding in a known sine-wave input at the profX'r 
frequency to give the same output jvading. The mean- 
squari' voltage per unit fretpiency was then obtained l>v 
dividing ly the effective band width. The eff’ectivt' band 
width had a valu(' of 5. 36 cps and was defined as tlie rc'c- 
t angular band width having the same area as the experi- 
nn'iilally d('t(‘rmiu('d band shaf)c. results present('d in 

figure's 6 and 7 are normaliz(‘d and 

and the' wave' number Ay is given by 


Over most of tlu' range tlu' accuracy was on the ord(*r of 
±10 perce'iit. The accuracy is somewhat less at tlu^ (wo 



Fi(;uK !0 5.- Distribution of turbulent she^ariiig stress and shear 
(M)rn*lation coefficient. 
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extrcmos because of large-amplitude fluctuations at low 
frequencies and because of low sigiial-to-noise ratio at the 
high frequencies. The error due to finite length of wire 
increases as the scale of the turbulent motions decreases and 
becomes significant at the higher wave numbers. However, 
because of the lesser accuracy of measureuKuit in this range 
no wire-length corrections wore applied. The overall value 


of serv(‘d as a check and was in good agreement with that 
cahnilated fi’om the spectrum. 

The same signals as those involved in the measurement of 
tlie overall turbulenc'e shear were passed through the same an- 
alyz(T to obtain the spectra of try shown in figure 8. Values 
below 10“^ should l)e taken with some reservation since they 
involve^ th(' small diff(‘jH‘nces of two nearly equal signals. 
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Figure 8. — Spectra of turbulent shear stress. 


The spectra of shown in figures 9 and 10 were calculated 
from the data taken during the measurements of the shear 
spectra, ddiis method, in principle, is perhaps not so accu- 
rate as a dire^ct measurement of the r-spectra, but a check on 
the reliability of the measurement was afforded by comparing 
the w-sp(H*tra calculated from the same data with that di- 
rectly measured, and in general the agreement was good. 

The (lata for the shear and spectra arc not normalized 
so that 



and 


PROBABILITY DISTRIBUTION, SKEWNESS, AND FLATTENING FACTOR OF 
u-FLUCTUATIONS 

The distribution of amplitude of ^^-fluctuations was de- 
termined by passing the voltage signal through a conven- 
tional gate circuit which conducted when the signal was 
above the level of the gate. The signal was preamplified so 
that the gate width could be considered neglible with respect 
to the root-mean-square value of the signal and the pulses 
produced by the gate were a series of square waves of con- 
stant amplitude of varying duration dependent on the time 
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Figuuk 10 “(Comparison of measured si)ectrum with that calculated from spectrum using isotropic relation; y/5^0.58. 


varied with the hn-el of the gate and ranged from 1 minute 
at the r)0-percent position to 5 minutes at the edge where the 
pulses were short and infrequent. 

The midpoint was d('fined as the position about which the 
av(T*ag(‘ signal was zero and is obtained from the first moment 
of the probability distribution, which by definition is 


^ 00 

= I uP(u) du=0 

v/ — 00 


The distributions of the probability chuisity at various 
positions across the boundary layer are shown in figure 11 
and are expressed in values of u relative to ?//, so that 


At most of the positions the distributions closed in rather 
well and permitted the calculation with reasonable accuracy 
of the skewness and flattening factors, where 


c* 

O f- 


j u^P(u) du ^ 
Q'” u^P{u) rfwj'' 


j u^Piu) du — 


At y/6—O.S and ?//5=1.0 the limbs of the curves were uncer- 
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tain, and no attempt was made to calculate the higher 
moments at these positions. Tlu‘ flattening factoi’ was also 
measured directly. Sinc(^ this involves the nondimeiisional 
ratio of the same voltage signal, that is, 

it was not ne(*essary to calibrate' the liot-wire, hut the cor- 
rection for the difference in gain involv(*d in the i)i‘ocessing 
of the divided signals had to he determined hy calibrating 
with a known sine-wave input. Tlie fourth power was 
measiu'cd l)y feeding the signal from the compensated ampli- 
fier into an instantaneous squaring circuit, which consisted 
of a series of diodes properly biased to give the sc|uare of the 
input (ref. 12), and then squaring again with a thermocouple. 
The distribution of flattening factors across the houndary 
layer measured in this manner is given in figure 12. 
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Figure 11. — Distributions of probability density. 
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Figure 12. — Distribution of flattening factor. 


It is seen from figure 11 that from yj6=0A to the wall the 
usual safety factor of three times the loot-rnean-square valuii 
for determining the ojierating point of the equipimuit is quite 
satisfactory. In the outer region of the layer, at ?//6=().8, 
this factor is exci'cded in the direction of nc'gative values of ?/. 
Here m(‘asur('inents were taken at varying gain-('ontrol 
settings to make sure tliat tlu' equipnu'nt was not 
ov('rloading. 

measliremp:nt of derivatives 

1'he usi'fuhu'ss of tlie concejit of local isoti’opy li(‘s in tin* 
ho[)e that the rnXv of turbulent-cuiergy dissipation in shear 
flow can b(^ given in terms of tlu‘ dissi])ation (expression for 
isotro()ic tinhuliMice. This would simplify tlu' experimental 
])ro(‘('dm“(' f)V ii'quiring only thi* measurement of the mean- 
square derivative' of the ?/-fluctuation with respect to j. 
A fuller discussion as to tlu' adaptability of the concept of 
local isotrojiy is given in the S(‘(‘tion entitled ^‘Energy 
Balance/' but it may be stated that the use of the isotropic 
ndation in eh'tc'rmining the dissijiation is inadequate', 
esp('cially in tlu' le'gion near thi‘ wall. It thus Ik'couk's 
necessary to iiK'asure all of the following nine mean-square' 
d('rivativesap|)('aring in the dissipation te'i'm of tlu' turbulence 
energy e(j nation 

(i)+iy+(fy+(fy 

In the present investigation, because of the practical 
limitations of hot-wire techniques it was ])ossible to measure 
ordy the first five of the above terms. The first three we're 
evaluated by taking the time derivative of the signal by 
means of a resistance-capacitance network and obtaining the 
mean-square value with a thermocouple and meter. The 
practical us(' of a differentiating circuit involves a comiiroinise 
betwf'en attenuation and the extent of linear response' with 
frequency. In the present measurements the differentiation 
was linear to 10,000 c‘ps and was down 12 percent fjoin 
linearity at lfi,()()0 (‘ps. This was considered adequate for 
the frequeiK'h's under invc'stigation since no significant change 
in tlie measun'd values was obtained by varying the cut-off 
filt(‘rs in the amplifier. The tiiiu* derivative was convi'rled 
to the sjiace dc'rivative by assuming the ac'curacy of the 
space-t imt' transformation 



The remaining two terms were obtained by measuring the 
correlation coefficients Uy and for a small distance of y 
and z and using the relations given by Taylor (ref. 17) 

,, u{y)u{y+y')_^ _ 1 (y'V 

u'{y)u'(y+y'} 2u^\^y/ 

y'-^O 

and 
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«'(2K(0+0') 

s'-^O 

The conelation coofficionts were measured by the method 
deseribod in reference 18 for varying distances of ?/' and 
ranging from 0.004 to 0.04 incli. 

distributions of the various derivatives a(*ross the 
boundary layer are give!) in figures 18 and 14. 

An at tiunpt to assess the validity of the space-time tratis- 
formation was made at y/5 = 0.05. In figure 15, the longi- 
tudinal correlation (coefficient Hx obtained by measuring the 
correlation lu‘ tween valiu's of u at the same instant for 
varying distan(“es in the x-direction is compared with that 
calcu!at(‘d from the speedrum using the Fourier transform 



Figurk 13.— Distribution of dissipation derivatives near wall. 



8 

Figure 14. — Distribution of dissipation derivatives away from w'all. 


7?i= J (A'l) cos {kiF) dki 

The latter gives the spatial correlation fi om a time con*ela- 
tion using the space-time relation. The microscale \x calcu- 
lated from the speed rum is dc^fimal by 

~,=l rir,^FJh) ilF 

Ar" -Jo 

It is seen that for values of x' less than 0.2 inch tlie two are 
in good agr(H‘ment. They yield tlu‘ saiiu' microscah*, and the 
space-time transfonnatioii is apparently valid for I lie small- 
scale motions n'sponsible for the dissipation. For lai'ge 
values of the two begin to diverge and the space-time 
relation becomes progressively worse for tht^ larger scale 
motions. At j:' = 0.4 incli the major (‘ontribution to the 
comdation comes from those wave numb(u’s below 1.3 per 
centimeter. It seems therefore that the ad(‘qua<*v of the 
space-time transformation depends on wave number and 
gives rise to the interesting speculation that tlie large-scale 
motions have tlieir own characderistic velocities different 
from the mean speed. Still (doser to the wall, the space-time 
assumption does not become a serious obstaede to the calcula- 



Figure 15. — Comparison of longitudinal correlation (‘oeflicient 
obtained directly with that calculated' from spectrum of ((2; y/5 = 0.05. 
Curve for X, also calculated from spectrum. 






12 


REPORT 1247 — ^NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


lion of t he dissipatioii since the magnitiRle is ^iveii mainly by 
derivatives in tlie ?/- and , e-directions. 

Although no wire-length corrections were applied to any of 
the data, it should bo mentioned that an estimate of the 

correction was made for at y/6r^0. ()()'> l>y the method 

given in reference 19. It was approximately 10 percent and 
conseqiHuitly wire-length corrections may (‘onsidered 
negligil)le across most of the boundary layer. 

DISCUSSION 

INTERMITTENCY^ 

In shear flows that have a free boundary, it has been re- 
peatedly obstu’ved tl)at as th(‘ free strtuim is approaclual the 
turl)ulence becomes intermittent, that is, that for only a 
fraction 7 of the time is the flow turbulent. This on-and-off 
chara{‘ter of the turbulcmce lias bi‘en (hdinitely establislied as 
being a manifestation of the irr egular outline of the boundaiw 
layer- as it moves downstr‘eam. The inter-mittiuicy is easily 
observed hy oscilloscope recoi'ds of the i/-fl net nation in the 
outer r-(‘gion of the Imundaiy layer-, and the recor'ds can be 
used to give a quantitative estimate of th(‘ factor y and to 
discei'ii some qualitative asp(M*ls of the flow, Kepresentative 
sections of oscilloscope r*ecor-ds taken at various positions 
across tbe boundary layer are given in figure 16. It is seen 
that in the outer region of the layer ?//5>0.4 ther-e are inter*- 
vals of time when the flow is not turbulent and that this time 
increases with inci’easing distance from the wall. Thus, the 
outer regime is divided into a turbulent part and a relatively 
non turbulent free-stream part, and the liot-wire at a given 
position resj)onds to alternate turbulent and non turbulent 
flows as the patti^rn is swept downstream. A vorticity meter 
(i‘cf. 20) which resi)onds to the vorticity that exists only in 
the turbulent region and suppr-esses tin* low-frequency 
fluctuation was used to obtain a sharper division between the 
turbulent and non turbulent regions. T^qjical samples of the 
vorticity trace o)x at different positions across the boundary 
layer are shown in figure 17. 

* Mrjusuroriu'nt.s by Corrsin ami Klstlor (n'f. 21) show tlu;t t)K‘ iiitcrinittcnl r<‘p:ioii in 
a boundary lay<>r ovi'r a roush wall ha.s tlw same general Ixdiavior. 


.4 V'*'' 










.86 



Th(‘ fir-st measurements of the intei*mitt(‘ncy factor y woi*e 
those obtained by Townsend (iTf. 3), where y is taken as the 
ratio of the flattening factor (of u or du/dt) in the wholly 
turbulent region to that in the iutermittcuit j-egion. Fr’om 
figure 12, it is seen that the flattening fardors in the i*(*gion 
near the wall ///6<C0.4 are quite constant with a valiK' cor- 
i-esponding closely to the Ciaussian value of d.O, By con- 
sid('ring the iiU tu’initteiicy as an on-off proc(‘ss tlu' value of 
7 is given by 


Th(‘ int(‘rmittencies calculated in this manner together witli 
those calculated from the oscilloscope n^airds of the -fluc- 
tuation and the vorticity are given in figure 18. Values of 
7 obtain(‘d from the flattening factor are (onsistent with 
those obtaiiU‘d ft‘om the film only up to y/d- OA). As sinui 
in figure' 12, tlu‘ fiatt('ning factors reach a maximum and tlum 
begin to di'cnaisi*. This is not too surprising bec-ausc' for 
low valiu's of 7 tile turbulence of the fi-('<‘ str(‘am would 1)0 
expected to make itself felt. The result is [U’obably a weight- 
ing of the jirobabilitv density for tlu' turbulence of the fn'(' 
sti-eaiii with tliat within the boundary, which dc'pends on 
the ri'lativi' tui-buhuua^ levels and the <h‘gr('(' of inti'rmittc'ncy. 
The curve' in figui-e* 18, which clos(4y nqire-seuits tlu' variation 
of 7 wuth 7//5, is a Gaussian integral curve given by 

7 = 7 ) (1— erf 0 

w hei'(‘ 
and 



44 le standard dc'viation a- is 

Such a distribution indicates that the instantaiu'ous j)osi- 
tion of the edge of tlu' boundary layer lias a random cliaracter 



Fkiuke 10.- -Turl)uU‘iice velocity u in boundary layer. Ui, 50 feet 
})cr second; timing dots, GO per second. 


Figukjo 17. — Turbulence velocity w* in Ixjundary lay(‘r. I\, 50 f(‘ct 
per second; timing dots, 00 per second. 
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with a moaii position at y/d~0.7S. The edge rarely extends 
outside the region y/8^0A to 1.2. While the position of 
the edge fluetuates over a large fraction of 5, the mean veloc- 
ity is iu‘ar that of the free stream, l)eing down by at most 
about 15 percent. 

Several interesting feature's can l)e gleaned from careful 
study of th(' oscilloscop(^ records. Although it may lu' difli- 
cult to sec^ from the short sections given in figure lb, it was 
noticed that in the strongly intermittent region, for example 
at ?//5=^0.8, the trace had somewhat of a square-wave apj)ear- 
ance. Increasing vi'locities are in the direction of the timiiig 
signal. The non turbulent regions set'in to be at a constant 
level corresponding to that of tlu' fn'c stivam, while the 
turbulent regions are seen to be centen'd about some lower 
level. The difference tx'tween the veIo(*ity of the outside 
potential flow and that existing in the turbulent retrion 
seemed to depend on how far past the measuring position 
(h(‘ instantaiu'ous edg(‘ of the layer extended at the particular 
instant. At ?//5=-1.0 and 1.2 there is very little evidence of 
the shift b(*caus(^ the edge dot's not exti'nd to any great dis- 
tanc(' Ix'yond these positions. Th(*re is then a step as well 
as th(^ on-off process which was a.ssumed in the calculation 
of 7 from the flattening fa(*tor. '^Du' agrt'ement with the 
valiK's obtained from tlu' film is {irobably due to the lack of 
sensitivity of tlu' flattening factor to this shift. However, 
the sk('wiH'ss factor can be exp(‘(*ted to be extr('mely s('nsitiv(*. 

A sket(‘h of tlu' boundary la\ (‘r is givt'u in figure 19. The 
boundary betw('en tiu' turbulent flow and the free stream is 
quiti' sliar{), and the ])roperti(‘s of tlu* slu'ar layer are com- 
paratively distinct from thos(* in the free stream. The 
boundary layer travi'ls downstream with an outline (*on- 
stantly changing in an irregular manner, and the intermit- 
U'ncy is (‘haracU'rized })v a large-scale diffusion firocess, 
carrying with it small-scale turbulent motions. From the 
film taken at ;y/5^0.8 it was also noti(*ed that the average 
frequeru'y of the occurrenci* of periods of no turbulence 
seenu'd to be ap[)roximat(‘Iy 100 per second, and, since the 





FicaiKK JS. — Distril)ulion of interniitU'ncy factor. 


pattern is moving with nearly the free-stream speed of 50 
feet per second, a rough estimate of the average wave length 
for the irregular outline of the boundary layer would appear 
to he appi'oxiniately 25 or 14a-. 

The effect of the intermit tency on tlie probability density 
is clearly s(‘en from figure 11. In the nonintermittent 
region tlu' distributions are very nearly Gaussian and the 
values of skewness and flattening factor calculated from the 
measured |)robabilit ies are given in the following table: 


y/5 

F S 

0.001 

2. :)2 0. 00 

.Of) 

2.1f^ -.OM 

.2 

2.02 -.OS 

.4 

3.19 -.20 j 


Calculations of the higher moments of the prol)ability 
densities from the distributions tend to be inaccurate because 
of tlu'ir emphasis on the higher values of ?/. However, the 
values of flattening factor agree fairly well with the directly 
measur('d vahu's. Values of skewness slmuld be taken with 
some reservation because they are ('xtremely sensitive to the 
accui'acy of tlu' midpoint. Also in the region m*ar the wall 
where turbulent fluctuations are large the nonlinear responses 
of the liot-wir<' tends to skew the signal. At 7//5 = 0.8 the 
probability d(*nsity is very strongly lU'gatively skewed 
because of the lower velocities within the turbulent regions. 
At the same time the large percent of time that the nontur- 
bulent regions exist causes the maximum to be displaced to 
the posit ivc side of the midpoint. The distribution at y/8 =-1.0 
has a flattening in qualitative agreement with the trend of the 
flatti'iiing factor to rc'ach a maximum and then to decrease. 
The distributions of skewness show the same general trend. 
The maximum of the probability distribution for y/8=l.O is 
displaced slightly to the left of and it is uncertain 

whi'ther this is expc'rinu'iital or a conscqiu'nce of the weighting 
effect previously mentioned in connection with the flattening 
factor. 

By assuming that the free-stream regions contribute little 
to the measured mean-turbulence quantities, an allowance 
may be made for the effect of intermittency by dividing by 
the factor 7. The distributions of turbulent energy and 
shear stress divided by 7 are given in figure 20. Tlie dis- 
tribution of turbulent energy within the bounded flow is 
strikingly similar to that for channel and pipe flow (refs. 4 
and 5) as is the distribution of turbulent shear stress which 
is approximately linear from ?//5=0.1 to ?//5= 1.0. 

Intermittency is, of course, absent from fully developed 


Boundary between layer 



Won 

Fkjukk 19. Sketch of boundary layer. 
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turbulent flow in pipes and channels because there is no free 
stream. It is well-known that the mean velocity for pipe 
and channel does not deviate as much from tlic logarithmic 
distribution as does that of the boundary layer, and a major 
factor in this discrepancy is the intermiUency. The in- 
fluence of the intermittency on the mean velocity distribu- 
tion is difficult to ascertain, being in the nature of a com- 
plex time-averaging problem. Although something may be 
doiie in a semiempirical fashion by introducing 7 into the 
mixing-length theories, it contrii)utes little to the basic 
understanding of the problem. In tigure 21 the measured 
mean velocity is plotted in the form sugg(*st(Hl by the veloc- 
ity-defect law, and it is seen that tlu^ chanu'teristic loga- 
rithmic, law exists only in a limited range fiom about y 5^ 
0.01 to yjd^O.2. The direct viscous dissi])ation (fig. 22) is 
negligible at ?//5 = 0.01. This corresponds to a value of 
27. G, which is in good agreement with the usually observed 
value of rU) for l>ipe and channel. An Intercast ing ])oint to 



I intTRE 20. — Turl)uknt shear stress and Uirt ulent energy divided i>\ y. 



Fir.uRK 21. — (kunparison of measured mean velocity distribution with 
velocity-defect law. 


note is that in the range of ?//5 — 0.1 to t//6=0.8 the ratio of 
shear stress to turbulent energy is approximately constant, 
as is also e\ idenced by the constancy of the shear correla- 
tion coefficient in figure 5. This would tend to indicate 
some degree of ‘dvarman similarity” existing in the range 
where the distribution deviates from the logarithmic law. 
In addition, Kai nnin similarity reciuires that the large-scale 
motions responsifde for the sliear stress be fret* from tlie ef- 
f(‘cts of viscosity, whicli, as will be seen in the seclion on 
spectra, is mu(“h more justified in the region where the 
logarithmic law is not obeyed than in the r<*gion wh(*re it is 
obeyed. Although any correction for the effect of iiiU'rmil- 
tency is in tlu* direction of minimizing this deviation, it is 
pointed out that the mean-velocity distributions for pipe 
and (‘haniu*! also deviate in this region. It is difficuilt to 
atta(‘ii any particular significance to the logarithmic law 
which exists for so limited a range of the boundary layer, 
(‘xcept that it is a J’cgion where the direct influence of the 
wall may still he present. It is not apparent that a degree 
of validity can be assigned to it on the basis of the various 
forms of tlu* mixing-length theory. With respect to the 
juean velocity, tlu* boundary layer can be divided into three 
regions, namely, a viscous i*egion exttmding to a value of 
= an inlermediatc region wdiere the influence of the* 
wall still exists, and an outer region characteriz(*d by 
intermittency. 



S 

Figure 22. Comparison of direct viscous dissipation with production 
of turbulent energy near wall. 
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ENERGY BALANCE 

After the usual boundary-layer approximations are made, 
the equation expressing the energy balance for tlie turbulence 
at a given cross section in a two-dimensional boundary layer 
is given by 




Tho resi)ective terms from left to right have the following 
physical interpretation : 

(1) Production of turbulent energy from the mean motion 

(2) Turbulent energy diffusion 

(3) Pressure diffusion 

(4) Convection of turbulent eriergy by the x-coinponent of 
the mean motion 

(5) Convection of turbulent energy by the y-cornponent of 
the mean motion 

(6) Dissipation of turbulent energy 

As discussed by Corrsin in reference 22 the viscous term 
in the energy-balance equation is not in a strict sense the 
true dissipation term. However, since for most of the 
boundary layer except for the region close to tlie wall, as well 
as in the integral across the boundary layer, it reflects the 
true dissipation, it is herein referred to as the dissipation 
term. 


convection of turbulent energy is negligible near tlie 
wall and is significant only in the outer region of the layer 
wheie similarity on the scale of 8 lias been experimentally 
observed, (Consequently, by using the continuity relation 


and 


da; ^ dy 


6oca-o'/-^ 


cept the diffusion terms which arc treated as one and obtained 
by balancing the equation. The distributions of the various 
energy terms ai'e given in figures 23, 24, and 25 in nondimen- 
sional form using 8 and Lq. The turbulent shear stress in 
the region 0<^y/5<^0.05 was obtained by subtracting the 
viscous sliear stress (calculated from tlie mean-velocity 
gradient) from the constant value of total shear stress as 
given by tlie sliear at the wall. The second derivative term 
is important only in the region near the wall. Its calculation 
is quite uncertain and is given mainly to show its order of 
magnitude, finite intercept at the wall is a consequence 

of the arbitrary linear extrapolation of the turbulence inten- 
sity to zero. The dissipation term is the most difficult to 
determine because of the importance of the region near the 
wall. I he applif’a ( ion of tlie concept of lo(’al isotropy is 
strongly conditioned by the local similarity of the eddies 
responsible for the dissipation, the verification of wdiich is 
dependent on the dissipation derivatives obeying the iso- 
tiopic relation: 



I’HiUHK 2.^. l^rodiu'l ion and dissipat ion in n'gion near wall. 


the convection terms can bo transformed into the single term 
which may be written in the nondimensional form 



wdierc is t lie distance from the virtual origin of the bound- 
ary layer. The dissipation term is not in a form tliat lends 
itself to measuriunent by hot-wire techniques, How^ever, it 
can be rewuatten in tlie following more suitable form: 


All of the pertinent term's in the energy balance can be 
calculated from the ineasurcinonts previously discussed ex- 



tuajiiio 24.' Production and dissipation in region away from wall 
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Kroni fijiurcs 115 jitkI 14 it is s(mmi lluit only in tin* oult^r 
rofiiou of the boumlary layer, ?//5>0.7, is this eonditioii satis- 
fied and that the derivatives become increasinydy divergent 


/du\^ , /duV 

as the wall is approached, d he terms \c)y) U 2 / 

equal across the layer, and at ?//5=0.00.) th('ir ratio to 


dx) 


is 10 times that pveii by isotropy, ddiis illustrates tlie very 
small scale naturi* of the turbulent motions in the ti*aiisverse 
dire(*tioiis as com])ared with the* loiij^itudiiial in the region 
close to the wall. The inade(iimcy of local isotropy is strik- 
ingly seiui fi’om the overall ener^j^y balan(*(‘ ol)taiiied by inte- 
grating the eiiergv (‘(piation across the boundarv layer, live 
integrals of the diffusion terms and the second-derivative 
t(‘rm ar(^ (‘a(*h zero, and tln‘ total jii’oduction is balanced by' 
the total dissipation plus the total convection. Using the 

isotropic relation 

turl)ul(uit energy, the I’csult is out of balance, witb the dissi- 
pation being too low bv a. factor of 2.1. Heuice, in order to 
obtain an accurate measure of the dissipation, it is necc'ssaiy 
to measure all nine derivatives. Sin<a‘ it \vas possible to 
measure only tlu' first five, the riunaining lour were obt aim'd 
by the arbitrarv assumption that derivatives with respc'ct to 
y and are given by' the* isotroi)ic n'lations 



This considerably improved tlu' total I'liergy balance. As 
seen from figure 23, the dissipation calculated in this manner 
was obtained as close to the wall as O.OOo, and if it is 



FiccKh. 25. Diffusion and ('{UiviHUion. 


assumed in order to close in the dissipation curve that the 
dissipation is equal to the production in the region 0<?y/5< 
0.01, the disparity is 15 percent. This can represent a sig- 
nificant diffenuice in the point-to-point balance as manifested 
by th(' diffusion term. When th(‘ diffusion term is compared 
with that obtained bv Townsend (ref. 0), the same general 
trend is observed although there is a considerable difference 
in magnitude. The necessary balance cannot be achieved 
bv the gain in energy due to diffusion in the region v//5>0.t) 
unless it is extrapolated to an unri'asonably long distance. 
Conseqin'Utly the conclusion is that the turbulent-energy 
(lissi[)at ion is greater than the prodiudion in a \ery^ thin 
i(‘gion iH'xt to the wall in order to provkh' the dissipation 
lU'eded. d'his means a diffusion of ('uergv toward tin' wall, 
d^ownsend s values for the kiiu'tic-i'iu'rgy'-diffusion term 
require a movenu'iit of kc'iu'th* (‘Ui'rgv’ away from the wall. 
Idierefore, if (uiergy goes toward the wall it must be due to 
tlu' action of pressure lorci's (pressure diffusion). ibis 
intensi' dissipative region for the turbulence appaiently 
coincides with the viscous n'gion for the mean flow'. I lie 
magnitude of the diffusion ti'rm illustrates the weakness of 
the mixing-length tlu'ory^ which has the implicit assumption 
that what is lo(*ally produced is locally' dissipated. 

^rhis flov' of ('uergy' against the energy^ gradient is ratlu'i* 
sur[)rising and contrary^ to what is intuitively ('xpected. 
It dot's not si'ein possible that such a result could Ix' brought 
about incorreet ly' b\' the questionable assumption math, 
regarding those dissipation derivatives that could not hci 
measured. Nt'ither does it st'cm reasonahh* that lln'St' 
derivatives can be distributed in such a manner as to makt' 
th(' production equal to the dissipation at all points across 
the layer. However, it is felt that the situation is compli- 
cated by tlu' phenomenon of intermittency and it would Ix' 
hel])ful to have a fuller understanding of the role that it 
plays. 

The ('iKMgy halauei' lor tlu' mean motion is conveiiic'iitly 
expn*ss(‘d in ti'rms of tlx' energy thickness 5, such that 



( 1 ) 


Similarity of tiu' mean velocity distribution ])crmits treating 
the integral in the expression for as a constant. If 5 is 
given bv 


then 


^Idc 




-l/5w, -1/5 
Jo 




K([tuUion (1) states that part of the loss of kinetic energy 
of the mean motion goes ilirectly into heat through the 
action of viscosity and the remainder goes into the production 
of turhidcnt energy. The energy balance was found to be 
satisfied to within a few percent and served as a welcome 
check on (he accurai’y of the measurements. The production 
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of turbulent energy ami the viscous dissipation are compared 
for Uic region near the wall in figure 22. This emphasizes 
the importance of the region near the wall. In fact, almost 
40 percent of the loss in kinetic energy of the motion is 
di^^ctly dissipated by viscosity in the region 0<y*<30, and 
of the remaining 60 percent which is converted into turbulent 
energy .'50 percent is produced in the same region. Thus, if 
the conclusion drawn from the turbulent-energy balance is 
correct, percent of the total dissipation (viscous and 
turbulent) takes place in y*<,'50. 

ENERGY AND SHEAR-STRESS SPECTRA 

J( is possible to ^oiin some further insight into shear How 
by examining the speetral distributions of turbulent energy 
and sh(‘ar stress. A significant advance in dealing with the 
(‘ri(‘rgy sjiectrum has l)een made in the domain of homo- 
geneous and isotropic turbuhuice (refs. 9 and 10). The 
basic concept underlying this advance^ is tliat emu-gy cuiters 
(h(‘ sp<‘(‘trum througli the large eddies and is then trans- 
ferred through th(‘ spectrum to the smaller eddies where it 
is finally dissipated. If the lower wave numbers are excluded 
ther<‘ exists a rang(^ in wdiich the eddies ar(» in a state of 
e(|uilibrium, governed by the rate at whi(*li they transfer 
end dissipate (Uiergy. When tlu^ Reynolds number is high 
eriough, inertial forces wdll predominate in the lower wave 
numbers of this equilibrium range, and a relatively pure 
Iranshu- r<*gion will exist. By dimensional re^asoning it can 
l)e shown that tlie spectrum will vary as in this range. 
Heisenberg (ref. 11) extended this (oncept by assuming 
that th(^ transfer of energy at wnve number k was caused 
by a turbulence friction produced by eddies wuth wnve 
number greater than He representerl the energy balance 
in the eejuilibrium rang(^ for homogeneous and isotropic 
t iirbidence as 



2{kJKit) dk' 


( 2 ) 


I Ih> se(*ond term wuthin the brackets represents a turbulent 
viscosity and X is an absolute constant. The solution of 
this (‘(juation gives the spectrum in the following form 

<3, 

wbcfc k, is a wave mimber in the intcrnK'diate range given by 





For low wav(> numbers E{k) varies as and for tbe high- 
wave-iumiber end where viscous forces predominate, as k~’’. 

Although equation (2) involves an assumed mechanism 
that may not entirely represent the facts, it seems to be a 
reasonably good approximation. Wher(! the Reynolds 
number is sufficiently high, there is evidence of a range 
and a transition to higher negative powers approaching —7 
as k increases. There is some doubt as to whether —7 is 
the corre(d, value in the limit. The concept embodied in 
equation (2) provides a rational liasis for an approach to 


the problem of shear turbulence. However, the extension 
to shear flow is complicated by the presence of such factors 
as production, diffusion, convection, and the absence of 
isotropy and homogeneity. Any conclusions as to the 
effects of the diffusion and convection are difficult to draw. 
It may be assumed that such effects are confim'd to the 
very low wave numbers which lie outside the (Hiuilibrium 
range. In addition, the convection may be considered 
negligible across most of the layer. An attemj)! to assess 
the influence of the production term in the e(|uilibrium 
range of the spectrum was made by fl’clien (ref. 23). By 
considering the influence of the mean velocity gradient the 
conclusion is rea«-hed that a range of /!•-' will ('xist in the 
wave-number region where k-^>^ normally exists when there 
is no gradient. 

An experimental test of theoretical predictions is rendered 
difficult by the fact that only th(^ one-dimi'iisional sjiectrum 
can be measured. In isotropic turbulence this is not too 
restrictive a factor because th(> transformation from the 
three-dimensional spectral function such as appears in 
equation (3) to the one-dimensional spectral function is 
known (ref. If). In shear flow this relation is not known, 
and one has to be content with the qualitat ive' infei i'iice that 
in some unknown manner the oiu'-dimensional siu-ctrum is 
still an integral effect of the thn'c-dirnensional. Despite the 
aforementioned complications the nu'asured sjiectra are of 
interest, and several interesting f(*atures can Ix' note'd. 

The spectra of at various cross-sectional jiositions an' 
given in figures 6 and 7. The trend, in going toward the 
surface, is for the higher wave numlH'rs to liai'e a gi'catt'i' 
percentage of the turbulent energy. This is in accoid with 
the trend of the shear-stress spectra (fig. 8) also to extc'iid to 
higher and higher wave numbers as the wall is af)])ioached. 
It is noted that nearly all of the turbident energy lies within 
the stress-producing range. The spectrum of slu'ar stress at 
»//5=0.2 shows an increase over that at ?//5=0.0,") in tin' lower 
wave numbei-s. This may be indicative of the influence of 
the wall becoming negligible at y/6=0.2. For 7//5>().2 the 
decrease in shear stress across the boundary la v('r takes j)lace 
for the entire spectral range, ^ill of tin* eneigy specTra 
indicate the e.xistence of a region varying as at the high- 
number end, but because of the lesser accuracy of the nu'as- 
urements in this range no direct comparison has bec'ii made. 
In the outer region of the layer at y/5=0..58 and 0.8 tlu'ic' is 
an extensive region where the spei'lrum of 7? varies as 
corresponding to the inertial subrange. The ('fleet of inter- 
inittency in this region is difficult to ascertain exce])t to say 
that it may be confined to the low-wave-numbc'r end of the 
spectrum. Also the maximum in the spectrum at 1 .0 is 
apparently a consequence of intermittency. 'rher(' is a 
gradual transition in the shape of the spectrum from y!8—().r,H 
to y 1 8 = 0.05 where there is a wav(>-number range with the 
slope ki~' as predicted by Tchen. The small-scale natun' of 
the turbulence near the wall is shown by the si)ectnim of 7? 
at i//5 = 0.0011. The dip in the spc'ctrum at llu' low-wave- 
number end as indicated by the dashed curve may he diu' to 
experimental error. However, it was rep('atahl(' and at ten- 
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lion is drawn to it bocaiiso it may bo a re^sult of sorno charac- 
toristic phoiiomonon asso(*iatod with tlio laminar siiblayor. 

Empirically, tho isotropic relation for tbe dissipation is in 
fair agreement with that given in figure 24 from yj6 = 0.i)5 to 
1.0. By assuming that the se(*ond moments of the spi'ctra 
(fig. 26) ai‘e a fair n>presentation of the total dissipation, it 
is seen that in going toward tbe wall tlu> shear penetrates 
deeper and deeper into the dissi])ation spectrum. For 
exam])le, at ?//5^^^0.r)8 tb<^ wave-numlier rang(^ up to when‘ 

fu.{k\) lias fallc'ii to ]()“■' coiitaiiir^ approximately 10 per- 

cent of the dissipation. At ?//5- 0.05 this figure is 80 percent. 
It is r(‘asonable to expect this ti(*nd to (*ontinue, with the 
i-esult that vmy close to the wall the turbulence produced is 
directly dissipated and their is no significant transfer of 
eniM'gy along tlu‘ spectrum. This conclusion is consistent 
with llu‘ small-scale turbulence ni^ar tlie wall and the rajiidly 
di'creasing speidrum at 7//5=^0.001 1 . 

The shear spectrum is a direct test of local isotropy, and 
it is evident that the range of wave numbers foi* which local 
isotropy exists becomes progressively smallm* as the wall is 
approached. Tlie transfer range, as evidenced by has 

not yet become locally isotropic since the shear spectra 
approach zero at some higher wave number. However, the 
eiu'rgy spectrum is ajiparently insensitive to tlie small amount 
of shear stress whicdi does exist. In figures 9 and 10 the 
immsured spiudra of V' at y/5-=0.05 and 0.58 are compared 
with those calculated from the measured spectra using the 
isotroi)ic relation 


m) 


A. 

2 2 dk, 


juik;) 


whi(‘h, strictly speaking, is valid only for the range of local 
isotropy. It is sihui that the measured spectra agree with 
the (‘alculated at the higher wave mimbei*s. This is con- 
sistent with th(‘ range' of lo(‘al isotropy indicated by the 
spectra of slu'ar, and th(' difference between the ovi'rall 

and V* is confiiual to tlu'- low(‘r wave numlx'rs. 

From the turbulent-energy ecpiations for the individual 
v(‘locity compoiu'nts it is s(‘en that llu' ti'i'ins responsible for 
th<^ i‘xclning<' of eiu'rgy among th(‘ various coniporu'uts are 

till' ])r('ssuj‘e t i'i'ins pi Jh ^^d(1 pi ihe speiural ranges 

of such transfer is unknown but a naisonabli! speculation 

would be that the i)ressure term pi res])onsibh> for the 

(raiisfi'r of eni'igy to v* I'xists in thi' same wave-numl>ei 
range as the shear-stress spe(*tra. There is little^ evidence 
in t he measured r sj)e(*trum of a range of as observu'd 

for the spectrum of u'. This may be a conseiiui'iu'e of the 
spectral distribution of the pressure term being different 
from that of tho shear stress with a weighting mori' to the 
higher wave number end of this range. 

CONCLUSIONS 

The following conclusions were obtained from an investi- 
gation of characteristics of turbulence in a boundaiy layer 
with zero pri'ssure gradient. Here y is the distance normal 
to the surface measured from the surface, y* is ?/ times the 
square root of the shearing stress at the wall divided by the 
density of air over the kinematic viscosity yUrjv (where 
f ^ is boundary-layer thickness. 




CHARACrJ’ERISTICS OF TURBULENCE IN A BOUNDARY LAYER WITH ZERO PRESSURE GRADIENT 
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1. Th(^ turbulent boundary lay<T (ran be (ronsidored to be 
dividcMl into three regions: A viscous region extending 
beyond the laminar sublayer to a value of = an inter- 
mediate region from ?/* = 30 to approximately y/6 = 0.2 where 
the influence of the wall still exists, and an outer region 
cliaracrterized by the phenomenon of interrnittency. 

2. The turbulent boundary layer has a sharp outline 
constantly changing in an irregular manner. The turbulent 
r(‘gions trav(‘l with a v(‘l()city low(‘r than that of tlu* fr(‘(r 
strc'am. A rough (‘stimaO* of the average wav(^ l(*nglh of this 
iiT(‘gular boundary would be aj)proximately 26. 

3. Th(‘ })osition of the edge of (he laycT with resj)(H*t to 
the surface is approximated closely by a Gaussian distribution 
centered at 0.786 with a standard deviation (Kpial to 0.146. 

4. ])robability density of the ^-fluctuation is v(u-y 
clos(* to Gaussian in tlu* region 7//6<^0,4. 

5. 'rhe imj)ortanc(^ of the iTgion near tlie vv^all has l>een 
d(*monstrated; in fa<‘t, a])proximately 85 p(‘rcent of the total 
dissipation (vis(*ousand turbulent) takes place witfiin y*<30. 
The ])roduction of turbulent (‘iiergy and t he t urbuhnit-energy 
dissipation leach a sharp maximum within this region, and 
th(‘r(‘ is an inward flow of emn’gy toward t lu* wall because of 
t lie action of j)ressure forces. 

0. Tlie coii(‘ci)t of lo(*al isotrop}^ is inadequate for obtaining 
the turbulent-energy dissipation, especially in the region near 
the wall. 

7. In the outer region of the boundary layer where the 
mean velocity gradient is small the existing theories for the 
eneigy sp<Hrtrum in isotropic turbulence can be applied. 
As th(‘ wall is approached the nonisotropy becomes significant 
witli a resultant change in the spe^ctrum. V(‘ry (‘lose to the 
wall, the turbulent energy produ(‘ed is directly dissipated with 
no significant transfer- of energy thmugh the spectrum. 

8. As closer to th(‘ wall as y/6 = ().05 it has benui shown that 
t h(‘ space‘-tiitie^ transfoemation is valid for the smaller eddie‘S 
r('S])onsi})le for tfi(‘ turbulent dissij)a(ion l)ut bea-omes invalid 
for the largen- (‘ddie‘s. 
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